Observation of Raman G-band splitting in top-doped few-layer graphene 
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An experimental study of Raman scattering in A''-layer graphene as a function of the top layer 
doping is reported. At high doping level, achieved by a CHF3 plasma treatment, we observe a 
splitting of the G band in the spectra of bilayer and 4-layer graphene (A'^ even), whereas the splitting 
is not visible in case of monolayer and trilayer graphene (A odd). The different behaviors are related 
to distinct electron-phonon interactions, which are affected by symmetry breaking and Fermi level 
position in different ways in the various A-layer graphenes. In trilayer graphene, a weakening of the 
electron-phonon coupling as a function of the Fermi energy induces a hardening of all zone-center 
in-plane optical phonon modes, like in monolayer graphene. On the other hand, in 4-layer graphene 
two distinct trends are observed in the G band as a function of doping, suggesting the presence of 
two different groups of electron-phonon interactions, like in bilayer graphene. 

PACS numbers: 



I. INTRODUCTION 

Since the demonstration of the isolation of a sin- 
gle atomic plane of graphite (graphene) on a standard 
Si02 1 Si substrate — , it has suddenly become possible to 
experimentally verify many theoretical predictions about 
the peculiar physical behavior of graphene-based sys- 
tems. Indeed, the linear E — k dispersion curves, and 
the consequent relativistic-like behavior of charge carri- 
ers in monolayer graphene, have been broadly confirmed 
by several experimental observations — i^. Interestingly, 
stacking a number N of graphene layers on top of each 
other can lead to new physical systems exhibiting com- 
pletely different properties. For instance, some peculiar 
gate-tunable electronic and optical properties have been 
recently reported in the bilayer case—"—. 

In fact, the interlayer coupling induces a gradual de- 
parture from the electronic bands of monolayer graphene 
— , until the bulk limit (graphite) is reached for N large 
enough. Therefore, there is a range of N where the phys- 
ical properties of graphene stacks are sensitive even to 
a variation AA^ — 1. For instance, a qualitative dif- 
ference between the transport properties of bilayer and 
trilayer graphene was evidenced in recent experiments — . 
Moreover, the optical absorption spectra are predicted 
to systematically vary with the layers number within the 
effective mass approximation for 1 < A^ < 6*, reflect- 
ing the gradual modification of the band structure. Here 
we show that a clear splitting of the Raman G band is 
observed for A^ = 2 and A^ = 4, when the multilayer 
graphene symmetry is broken by heavy doping of the top 
layer, whereas the splitting is not observed for A^ = 1 and 
A^ = 3. The presence or absence of splitting highlights 
different electron-phonon interactions, which are influ- 
enced by doping and symmetry breaking in distinct ways 
in the various A^-layer graphenes. The reported results 
confirm two experimental reports in literature about the 
G band splitting in bilayer graphene (obtained by gate 
field effect adding new important information such 

as a systematic study of the splitting as a function of the 



layers number and of the doping. 

Raman spectroscopy is a very powerful tool for study- 
ing graphene, yielding information on the electronic 
structure and on the electron-phonon coupling (EPC) in 
the material— ii^. In fact, this technique allows to clearly 
distinguish a monolayer and a bilayer from a few-layer 
graphene, through the analysis of the 2D band — , and 
to estimate the charge carrier density and type in mono- 
layer graphene from the spectral positions and relative 
intensities of the G and 2D bands — . The G band (at 
^ 1580 cm"^) is due to a first-order Raman scattering 
process involving zone-center in-plane optical phonons. 
In stacked graphene layers, the vibrations in different 
atomic planes can combine with each other in various 
ways, depending on the number of layers and on the 
stacking order. Moreover, the EPC is affected by the 
number and symmetry of the stacked layers, so that the 
G band may be used to study the effects of symmetry 
breaking and doping on the electronic and phononic prop- 
erties of multilayer graphenes. The knowledge of these 
effects is of fundamental importance for the development 
of graphene-based field-effect electronic devices. 



II. EXPERIMENTAL 

Graphene layers studied in this work were deposited 
on a 285 nm thick Si02 on Si substrates, by adhesive 
tape exfoliation of natural graphite. Then, the samples 
were analyzed by optical microscopy, in order to estimate 
the number of graphene layers composing the deposited 
thin flakes. It can be seen in Fig. [T] that the contrast 
(defined as 1 — Rq/Rs, where Rg and Rs are the re- 
flected light intensities from the Si02/Si substrate with 
and without graphene, respectively) measured on many 
semi-transparent fiakes increases in a stepwise manner. 
The analysis of the 2D Raman band confirmed that the 
lowest two steps correspond to monolayer and bilayer 
graphenes, indicating that few graphene layers can be 
counted by contrast analysis. Such a behavior is related 
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FIG. 2: Raman spectra of monolayer graphene after plasma 
treatments at various CHF3 flow rates 



13 cm ^) indicates that the unintentional doping level 
in the pristine samples was relatively low (about 1X10^^ 
cm-2)i8. 

Then, the samples were subjected to a radio frequency 
(RF) plasma treatment in CHF3 gas, and immediately 
(within a few minutes) re-characterized by Raman spec- 
troscopy. 



III. RESULTS AND DISCUSSION 



FIG. 1: (color online). Optical microscope image of one of 
the flakes analyzed in this work (top). The number of layers 
estimated by contrast analysis is indicated. In the bottom 
graph, the discrete behavior of the contrast as a function of 
the number of layers is shown. Experimental data, obtained 
at A=550nm, are compared to theoretical values. 

to the optical absorption of graphene, which was found 
to be directly proportional to the number of layers for 
N small enough Using appropriate filters in order 
to select the wavelengths at which the contrast varia- 
tion is high, we were able to distinguish up to 6 layers. 
Moreover, the experimental contrast values were checked 
by theoretical calculations within the Fresnel coefficients 
approach!^. After a preliminary annealing in vacuum 
(lxlO~^ mbar for two days) in order to remove possi- 
ble adsorbed impurities from the graphene surface, the 
first run of Raman measurements was carried out. Ra- 
man spectra were acquired by means of a Jobin-Yvon 
UIOOO Raman spectrometer equipped with a microscope 
(lOOX objective) and with an Ar-Kr laser, using the ex- 
citation wavelength A = 514.5 nm. The incident laser 
power focused on the sample was adjusted to be less than 
5 mW, to avoid any local heating effect. Various iV-layer 
graphene flakes were analyzed, displaying the standard 
G and 2D bands reported in literature. The D band at ^ 
1350 cm~^, related to lattice defects, was never observed 
in the experiments, confirming the good quality of our 
graphene samples. Moreover, the analysis of the mono- 
layer spectra (G peak at ~ 1582 cm'^ with FWHM 



A. CHF3 plasma 

Previous studies reported in literature— 1^ have shown 
that the following radical species can be found in a CHF3 
plasma: F atoms and CF^ {x = 1,2,3) radicals. In dry 
etching processes, F atoms normally act as the reactive 
species (responsible for the etching), whereas CF^ radi- 
cals are passivation precursors giving rise to polymer de- 
position. It has been observed that, in the case of CHF3 
plasma processes, the etch rate decreases with increasing 
flow (in contrast with what happens in the case of CF4 
plasma), due to the low concentration ratio [F]/[CFa;]— . 
This means that the action of F atoms becomes more 
effective with decreasing the flow rate, because the pas- 
sivating action of CF^ radicals is reduced. In our exper- 
iments, we have observed an increasing modification of 
the graphene Raman spectra, with decreasing the flow 
rate in the plasma treatment. 

In particular, we performed a preliminary study on 
monolayer graphene processed at various CHF3 flow 
rates. Fig. [5] shows that the plasma treatment induces 
a blue shift of both G and 2D peaks, which increases 
with decreasing the gas flow. Moreover, below a flow rate 
threshold (about 6 SCCM), two new peaks arise at about 
1350 cm^^ {D peak) and 1620 cm^^ {D' peak), which in- 
dicate the presence of defects in the sp^ G lattice. Also 
the effect on the Raman spectra of A^-layer graphenes, 
which is going to be discussed in detail in the following 
Section, was remarkably reduced with increasing the gas 
flow rate. 
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In analogy with the mechanism involved in the etching 
processes, the interaction of F atoms with the graphene 
surface is likely reduced at high flow rate, because of 
the passivating action of the CFx radicals. Therefore, 
the effect of CHF3 plasma on graphene Raman spectra 
can be ascribed to the adsorption of F atoms on the 
surface. With increasing the F coverage, the modifi- 
cation of the graphene properties becomes more impor- 
tant, passing from a p-type doping effect (blue-shifted G 
and 2D peaks) — to a structural modification {D and 
D' peaks) at a very low flow rate (less than 6 SCCM). 
The doping effect can be ascribed to an electron trans- 
fer from graphene to adsorbed F atoms, i.e. a mecha- 
nism analogous to the observed charge transfer between 
graphene and adsorbed K atoms At very low gas fiow, 
when the passivating action of the CF^ radicals is min- 
imized, chemical modification (fiuorination) of graphene 
may eventually occur, with a transition from sp^ to sp^ 
C hybridization similar to that observed in graphane for- 
mation by plasma hydrogenation— . This evolution may 
be analogous to the transition from semi-ionic to covalent 
C — F bonding observed in carbon nanotubes treated in 
CFi plasma—. 

Here, we discuss the results of processes carried out 
at a gas flow of 6 SCCM (pressure of 100 mTorr) for 5 
minutes, at RF power = 15 W. Such experimental con- 
ditions lead to a very high doping, without structural 
modification of the graphenes. Moreover, the symme- 
try of stacked graphenes is broken by the dipole moment 
generated by the charge transfer from graphene to the 
adsorbed F atoms. Therefore, the situation under study 
is very similar to that found in field-effect experiments, 
and the results here reported may be useful for the study 
of gated graphene-based devices. 



B. G band splitting and electron-phonon coupling 
in A'^-layer graphene 

We focus now on the effect of top doping on the Raman 
G band, which displays very distinct features depending 
on the number N of stacked graphene layers. The change 
of the G band induced by plasma treatment in the vari- 
ous cases is visible in Fig. [31 where the spectra at to and 
ti were acquired on the same substrate before and im- 
mediately after the treatment, respectively. The G peak 
of monolayer graphene is largely blue-shifted (to ~ 1590 
cm^^) and narrowed (FWHM ~ 6 cm^^). Both the ob- 
servations are consistent with an increase of the doping 
level, which induces a hardening of the mode, due to the 
non-adiabatic removal of a Kohn anomaly for zone center 
optical phonons , and a reduction of the linewidth, due 
to Pauli exclusion principle which inhibits phonon decay 
into electron-hole pairs when the Fermi level surpasses 
half the phonon energy ~. On the other hand, the bi- 
layer and 4-layer spectra display a very evident splitting 
of the G-mode, whereas the behavior of trilayer spectrum 
is similar to that of monolayer. These results have been 



confirmed on different flakes on the same sample and on 
different samples. 

Furthermore, we observed that the modification in- 
duced by the plasma treatment was not stable under 
ambient conditions, as the Raman spectra changed with 
the passing of time, slowly tending to their pristine form 
(Fig. [3]). Indeed, the initial condition can be restored by 
a vacuum annealing, so that the plasma treatment can be 
repeated for several times in a reproducible way. Such a 
reversible behavior is consistent with the absence of the 
D peak, which indicates the lack of structural modifi- 
cation of the material (in contrast with the case of the 
chemical modification obtained by plasma treatment at 
very low gas fiow). As previously discussed, the mono- 
layer spectra, showing the hardening of both the G-mode 
and the 2D mode, tell us that a strong p-type doping is 
achieved — upon plasma processing. Therefore, we were 
able to gradually vary the doping level on top of each 
A^-layer graphene, studying the effect for different values 
of TV. 

The effect of charged adsorbates has already been stud- 
ied on epitaxial bilayer graphene on SiC by ARPES mea- 
surements It was shown that the electronic bands 
of bilayer graphene are strongly affected by a potassium 
atoms coverage on the top surface, due to the n-type 
doping induced by the adsorbates. Indeed, because of 
the short screening length along the c-axis of stacked 
graphenes — , the majority of the doping charge is lo- 
calized in the top layer and a dipole moment is formed, 
giving rise to a symmetry breaking and to the conse- 
quent band gap opening — . In our case, an analogous 
effect is likely to occur: the doping charge density rapidly 
decreases going from the top layer to the bottom, thus 
altering the symmetry of the multilayers. The symme- 
try breaking can affect very much the physical properties 
of the system, as it eliminates the symmetry constraints 
which rule both phonon Raman activity and EPC. In 
these conditions, all phonon modes included in the G 
band can become Raman active and strongly mixed with 
each other—, and the EPC is affected by the modifica- 
tion of the electronic bands. 

In Fig. [31 the G band dependence on the doping is 
visible in all cases, but it looks different for different 
values. An estimation of the doping level may be ob- 
tained from monolayer graphene spectra, basing on the 
G peak position and linewidth, and on the intensity ratio 
of the G and 2D peaks but we have verified that dif- 
ferent monolayers can display slightly different G peak 
positions on the same sample, due to the difficulty of 
controlling the native doping of graphene in ambient at- 
mosphere-". Furthermore, monolayer graphene is likely 
to have a different reactivity with respect to that of few- 
layer graphene, as experimentally observed in the case of 
hydrogenation ^1^. Also in our experiments, we have 
seen that, after plasma treatment at very low flow rate, 
the D peak is hardly observed in the few-layers spec- 
tra, whereas the monolayer spectrum displays a very ev- 
ident D peak. Therefore, it seems not correct to extrap- 
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FIG. 3: (color online). Evolution of the Raman G band of A'^-layer graphene after CHF3 plasma treatment. The spectra at 
to were acquired before the treatment, while those at ti (1 < i < 6) were taken at various time intervals after the treatment 
{ti ~ 15mjn;t2 ~ 6/1; t-j, ~ 24/i; ti ~ 48/i; ts ~ 72/i; ta ~ 144/i ). The lines connecting peaks are guides for the eye. 



olate the doping values from the monolayer analysis to 
interpret the few-layer spectra. However, although a pre- 
cise quantitative estimate of carrier concentration in each 
case is not possible in our experiments, we can monitor 
the different behavior of A^-layer graphene Raman spec- 
tra with decreasing the top doping, starting from high 
doping levels (more than 1 x lO^^cm"^) as suggested by 
the analysis of monolayers spectra. 

The effect in the bilayer can be interpreted on the ba- 
sis of some recent literature. Indeed, the splitting of 
the Raman G band was recently observed in gated bi- 
layer graphene and ascribed to the inversion symme- 
try breaking induced by the gate field effect and to two 
distinct EPC involved in the G band. This is due to 
the fact that the G band of bilayer graphene includes 
two doubly degenerate modes, i?2g and £'„, which are 
symmetric and antisymmetric with respect to the inver- 
sion symmetry, respectively. Consequently, only the E2g 
mode is normally Raman active, unless inversion symme- 
try breaking switches the mode on, too. It has been 



shown by Ando and Koshino — that, in the presence 
of asymmetry in the potential of the two stacked layers, 
symmetric and antisymmetric modes are strongly mixed 
with each other and two peaks appear in the Raman spec- 
trum. Moreover, phonons can be considerably modified 
by resonant electronic interband transitions, when the 
asymmetry opens up a gap comparable to the phonon 
energy. Ab initio calculations have been performed to 
compute the G band of bilayer graphene in asymmetric 
conditions as a function of the carrier concentration in 
the top and bottom layer [ritop and ribot) predicting 
the behavior of the two modes which become Raman vis- 
ible at certain values of ritop and n{,ot- 

In Fig. 01 we report the analysis of the two Lorentzian 
peaks which can fit the bilayer G band at various stages of 
the plasma modification, for two different samples. The 
experimental evolution of the peaks position has been 
fitted by the theoretical curves obtained from refi^, as- 
suming that about the 85% of the total charge carriers is 
confined in the top layer, according to ref.--, and using 
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TABLE I: Irreducible representations of the IR and Raman 
active modes at the F point for A^'-layer graphene 



N 






Even 
Odd 


(iV - l)A2u (S{N- 

NA2 e ne' 


NAlg e NEg 

(N - 1)A[ ® ne' ® {N - 1)e" 



the total carrier concentration n as the free parameter 
in the fitting. Then, the other spectral features (inten- 
sity and linewidth) have been compared to the behavior 
predicted in refj^, considering the carrier concentration 
values obtained from the best fit of the peaks position 
curve. It can be seen that a good qualitative agreement 
is obtained for all the analyzed parameters. In particular, 
the threshold at about n — 1 x 10^^ cm^^, after which a 
steep variation of both the intensity and linewidth ratios 
is predicted by theory, is well reproduced by the experi- 
mental data, thus confirming the consistence of the car- 
rier concentration results obtained by fitting the peaks 
position curves. Therefore, the behavior of the bilayer 
G band is well interpreted in the framework of an asym- 
metric carrier distribution model. 

In order to interpret the spectra for > 2, it is worth 
to consider, first of all, the evolution of the G band with 
N as predicted by group theory. Indeed, applying the 
group theory to Bernal stacked graphene layers, the irre- 
ducible representations of the infrared and Raman active 
modes at the F point for A'^-layer graphene can be ob- 
tained ^'^1'^^ , as listed in Table H] The modes related to 
the G band are the Eg, Eu, E and E in-plane modes, 
while the A modes are related to out-of-plane phonons. 

Eg and E^ modes are found in case of inversion sym- 
metry of the system (N even) , whereas E and E modes 
appear for mirror symmetry {N odd). Importantly, only 
the inversion symmetry inhibits the Raman activity of 
antisymmetric modes. It can be obtained by ab initio 
calculations — that in trilayer graphene two E modes 
and one E mode can be found at the G-band frequen- 
cies, whereas in 4-layer graphene two Eg modes and two 
Eu modes vibrate at the G band frequencies. The in- 
version symmetry breaking, produced by the top dop- 
ing, switches on the previously Raman silent antisym- 
metric modes i?„ in N even-layer graphene, so that all 
the phonon modes become Raman active and mixed in 
A^-layer graphene for every N value. Moreover, like in 
the bilayer case previously discussed, the EPC can be 
strongly affected by the lack of symmetry constraints and 
by the change of the electronic band structures. 

The absence of splitting for trilayer graphene can be 
qualitatively interpreted considering the allowed elec- 
tronic interband transitions which give rise to the phonon 
energy renormalization (Kohn anomaly) when Ep 
(low doping). In Fig. [5] it is shown that both E and 
E phonons can couple with electronic transitions when 
the Fermi level is at the Dirac point, whereas in bilayer 
graphene only the symmetric Eg mode can efficiently cre- 
ate electron - hole pairs, due to energy conservation and 
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FIG. 4: Features of the two peaks fitting the bilayer G band as 
a function of doping. Circles and triangles are experimental 
data from two different samples, and lines are the theoretical 
findings of Ref.— . In (a), the peaks positions have been fitted 
by the theoretical curves, assuming that about the 85% of the 
total charge carriers is confined in the top layer. In (b) and 
(c), the width and intensity ratios of the low (L) frequency 
and high (H) frequency peak are reported. 



symmetry selection rules. Indeed, first principles calcula- 
tions 3& have shown that in trilayer graphene the phonon 
linewidths of symmetric and antisymmetric modes are all 
of the same order of magnitude, whereas in the bilayer 
the antisymmetric linewidth is two orders of magnitude 
smaller than that of the symmetric mode. Therefore, al- 
most the same EPC strength is expected for the three 
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FIG. 5: (color online). Electronic transitions allowed by sym- 
metry rules in bilayer (left) and trilayer (right) graphene when 
tF = O.The electronic bands are taken from Ref.*. The transi- 
tions indicated by dotted lines are not involved in the phonon 
energy renormalization, because their energy is much higher 
than the G band phonon energy (~ Q.19&eV), which is re- 
ported as a scale bar on the left. 



the band structure of the 4-layer, which can be approxi- 
mated by two bilayer type band structures — . In partic- 
ular, intraband transitions occurring between two sub- 
bands separated by /S.E ~ 0.2AeV, which corresponds 
to an experimentally observed absorption by IR spec- 
troscopy may give rise to an almost resonant coupling 
with phonons {Eph ^ 0.2eV), and to a consequently 
strong renormalization of phonon energy. This EPC 
is expected to have a maximum when the Fermi level 
reaches the high energy sub-band (see Fig. ^ , when the 
number of possible intraband transitions is maximized. 
Therefore, the doping value n,„, at which a minimum 
for the Raman peak position is observed, is likely to cor- 
respond to the Fermi energy touching the high energy 
sub-band at Ep ~ 0.2AeV. A theoretical analysis of the 
4-layer band structure and density of states in the pres- 
ence of asymmetric doping may confirm this hypothesis. 




FIG. 6: (left) Electronic intraband transitions in 4-layer 
graphene for ep ~ 0.24el/. The electronic bands are taken 
from Ref.~, and the G band phonon energy (~ 0.196el/) is 
reported as a scale bar. (right) Positions of the two peaks 
fitting the splitted G band of 4-layer graphene as a function 
of doping. The observed dip may be related to a strong EPC 
due to the intraband transitions shown on the left. 

phonon modes in trilayer graphene, giving rise to the 
phonon energy renormalization at Ep — and to the con- 
sequent hardening of all phonons with moving the Fermi 
level, like in monolayer graphene. This is a qualitative 
different case with respect to that of bilayer, where only 
the symmetric phonon energy is renormalized at Ep ~ 0. 
The presence or absence of the G band splitting in the 
two cases reflects the presence or absence of distinct EPC 
for phonons of distinct symmetry. 

The observation of splitting in the 4-layer spectrum 
suggests the presence of distinct EPC for the different 
phonon modes in analogy with the bi-layer behavior. In 
Fig. [B] we show the experimental position of the two 
lorentzian peaks fitting the 4-layer G band, as a func- 
tion of doping as obtained from the bilayer analysis. It 
is worth noting that a minimum can be clearly identified 
in the curve of the low frequency peak, indicating that a 
maximum is likely to occur in the EPC for some phonon 
modes at a given doping value (rim — 9X10^^cm~^). 
Such a feature can be interpreted by taking into account 



IV. CONCLUSIONS 

In summary, we have experimentally investigated the 
Raman G band for Ai"-layer graphene (1 < < 4) in the 
presence of high asymmetric doping, finding two different 
types of behavior. For N odd, the G band is always 
fitted by a single lorentzian peak, which is blue-shifted 
with increasing the doping level. This is due to a strong 
EPC for all phonon modes when Ep ^ 0, which decreases 
with increasing the Fermi energy. On the other hand, for 
N even, an evident splitting of the G band is observed, 
related to the presence of distinct EPC for phonons of 
distinct symmetry. In particular, in the 4-layer case a 
signature of the van Hove singularity at E ^ Q.2AeV is 
likely to be observed as a minimum of the low energy 
peak position. 

Insights into the electron - phonon interactions in iV- 
layer graphenes in the presence of top doping can be use- 
ful for the study of field effect graphene-based devices. 
Moreover, the CHF^ plasma treatment may be a pow- 
erful technique for the study of graphene in the presence 
of a coverage of highly electronegative atoms such as flu- 
orine. 

Finally, the variety of the results reported in the lit- 
erature about the Raman G band in heavily doped bi- 
layer graphene —liii^i^ suggests that the repartition of 
the additional charge carriers is not well understood in 
most experiments. Therefore, suspended graphene sam- 
ples may be a good testbed to further investigate the dis- 
tribution of the doping, especially in the bilayer case4Si^ 
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